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Abstract 

The coupled cluster method (CCM) is a well-known method of quantum 
many-body theory, and in this article we present an application of the CCM 
to the spin-half J1-J2 quantum spin model with nearest- and next-nearest- 
neighbour interactions on the linear chain and the square lattice. We present 
new results for ground-state expectation values of such quantities as the en- 
ergy and the sublattice magnetisation. The presence of critical points in the 
solution of the CCM equations, which are associated with phase transitions 
in the real system, is investigated. Completely distinct from the investigation 
of the critical points, we also make a link between the expansion coefficients 
of the ground-state wave function in terms of an Ising basis and the CCM 
ket-state correlation coefficients. We are thus able to present evidence of the 
breakdown, at a given value of of the Marshall-Peierls sign rule which is 
known to be satisfied at the pure Heisenberg point ( J2 = 0) on any bipartite 
lattice. For the square lattice, our best estimates of the points at which the 
sign rule breaks down and at which the phase transition from the antifer- 
romagnetic phase to the frustrated phase occurs are, respectively, given by 
^ ^ 0.26 and ^ f» 0.61. 
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I. INTRODUCTION 



Antiferromagnetic materials have conveniently been modelled, since the 
early work of Heitler and London, as a lattice of magnetic atoms upon which 
the active electrons are localised. Furthermore, the exchange interactions be- 
tween the electrons are conventionally described solely in terms of the spin 
degrees of freedom of the electrons. An archetypal model of this class remains 
the Heisenberg model in which only nearest-neighbour exchange interactions 
are included, and these are all taken to be equal. Although the Heisenberg 
model on the one-dimensional (ID) chain has been exactly solved many years 
ago by Bethe [1], it is still the case that relatively few other exact solutions 
have been found in the intervening 65 years or so to comparable models in 
higher dimensions or to models involving more complicated interactions, es- 
pecially those containing an element of frustration. 

On the other hand, various approximate numerical techniques have by now 
been applied to a large number of such magnetic lattice Hamiltonians. For 
example, many variational calculations have been undertaken, employing a 
wide variety of trial wave functions. Although these often give accurate upper 
bounds for the ground-state energy, for example, one often finds that the 
differences between the estimated energies for trial states of widely differing 
kinds are very small. Hence, predictions based on variational calculations for 
properties other than the energy, or to such questions as whether the exact 
ground state is ordered or disordered, are notoriously unreliable. 

As a common alternative one may perform exact diagonalisations on small 
finite clusters of spins drawn from the infinite lattice under consideration. 
However, even with modern computers, one considers clusters of A'^ spins with 
N < 36. Extrapolation to the infinite lattice, N ^ oo, then needs to be per- 
formed. While exact results from finite-size scaling theory are often of great 
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help in this regard, the extrapolation does need to be handled with great care. 
This is particularly true when using the results from finite clusters to make 
quantitative predictions for such quantities as the order parameter. Many 
wrong claims have been made in the past from an improper treatment of the 
very subtle phenomena which need to be taken into account in the extrapola- 
tions, as has been stressed and discussed with great care by Lhuillier and her 
co-workers [2]. Furthermore, one expects that such finite-cluster calculations 
will become less accurate the closer one approaches a quantum phase tran- 
sition between states of different quantum order, marked by a critical value 
of some coupling parameter, at which a correlation length characterising the 
order typically diverges. 

Results for much larger clusters are typically obtained by stochastic simu- 
lation of the many-body Schrodinger equation using various quantum Monte 
Carlo (QMC) algorithms. Where the basic spin-lattice Hamiltonian can be 
mapped onto an equivalent bosonic problem, as in the case of the Heisenberg 
model on a bipartite lattice, such QMC techniques can readily be applied 
to clusters containing several hundred or more spins, and very accurate re- 
sults thereby obtained. In these cases, such as the Heisenberg model on the 
two-dimensional (2D) square lattice [3,4], the QMC results can usually be 
considered as "exact" , with the resulting errors arising only, or largely, from 
statistical errors which are open to systematic reduction within the limits of 
the available computing power. 

What ultimately underpins these bosonic mappings, and what therefore 
makes such QMC simulations so readily attainable, is the knowledge that in 
some appropriate representation the multi-spin cluster coefficients describing 
the iV-body wave function are all positive-definite. For example, in the case 
of the Heisenberg model on a bipartite lattice, this information is provided 
by the well-known Marshall- Peierls sign rule theorem [5]. 
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Conversely, where such prior knowledge of the nodal structure of a many- 
body wave function is not exactly known, QMC calculations are beset by 
the notorious "minus-sign problem", and are usually then much less reliable 
or much more difficult to implement with known algorithms. A typical way 
that such complications arise in spin-lattice problems is from the introduction 
of frustration. This can arise either from the geometry of the lattice under 
consideration or from the introduction into the Hamiltonian of competing ex- 
change interactions. An example of the former is the basic Heisenberg model 
on a 2D triangular lattice; while an example of the latter arises from the intro- 
duction on a bipartite lattice of (antiferromagnetic) next-nearest-neighbour 
interactions in addition to the (antiferromagnetic) nearest-neighbour interac- 
tions of the pure Heisenberg model, resulting in the so-called J1-J2 model 
studied here. 

Relatively few QMC calculations on such frustrated spin-lattice systems 
have been performed. As a starting point they typically require a good trial 
wave function, in terms of which the true wave function is well approximated, 
especially for its nodal surface structure. In such calculations there can still 
be a considerable systematic uncertainty, beyond the unavoidable statistical 
errors, which arises from whether the simulations have eliminated the bias 
inherent in the starting function. A typical recent calculation of this type 
was the fixed-node Green function Monte Carlo method [6] simulation of the 
spin-half 2D triangular-lattice Heisenberg antiferromagnet by Boninsegni [7] . 
While undoubtedly representing a very ambitious calculation of its kind, the 
resulting prediction for the sublattice magnetisation, which is the simplest 
measure of the antiferromagnetic Neel long-range order in this system, seems 
to be clearly too high by comparison with the results from the best of the al- 
ternative techniques, including exact diagonalisations of small clusters [2] and 
the coupled cluster method [8]. Furthermore, even the resulting QMC esti- 
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mate for the ground-state energy of the triangular Heisenberg antiferromag- 
net gives an upper bound which is relatively easily bettered by the alternative 
techniques. 

For such frustrated and similarly "difficult" systems, predictions based on 
even very large-scale QMC simulations share, at least to some extent, the un- 
certainties discussed above for variational calculations. In order to overcome 
these uncertainties, therefore, there is a real need to apply alternative semi- 
analytical approaches, especially those that have the demonstrated power 
to provide accurate predictions for the quantum order and for the positions 
and critical properties of any quantum phase transitions. One such method, 
namely the coupled cluster method (CCM) [9-19], stands to the fore in this 
respect. It has long been acknowledged as providing one of the most powerful, 
most widely applicable, and numerically most accurate at attainable levels of 
computational implementation, of all available ab initio formulations of mi- 
croscopic quantum many-body theory. Furthermore, in recent years it has 
been widely applied to many spin- lattice Hamiltonians [8,20-33]. For exam- 
ple, very successful applications have by now been made to the solid phases 
of ^He [20]; the isotropic Heisenberg and anisotropic XX Z models in ID 
and on the 2D square lattice, both for spin-half systems [20-28] and higher- 
spin systems [29]; and the spin-one Heisenberg-biquadratic model on the ID 
chain [30]; as well as to such frustrated spin-half models as the J1-J2 model 
in ID (and 2D) [31-33] and the 2D triangular-lattice anisotropic Heisenberg 
antiferromagnet [8,26,28]. 

In the present paper we apply the CCM specifically to investigate phase 
transitions in the spin-half J1-J2 model on (the ID chain and, especially, on) 
the 2D square lattice. Our main aim is to use this model as an archetypal ex- 
ample for which no exact information is known for the nodal structure of the 
exact ground-state wave function, apart from the Marshall-Peierls sign-rule 
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results in the pure Heisenberg limit. There is particular interest in studying 
whether the sign rule is preserved when weak next-nearest-neighbour exchange 
interactions are included and, if so, whether there is a critical coupling beyond 
which the sign rule breaks down. We believe that the CCM is an excellent 
prima facie candidate for such studies, since in virtually all previous appli- 
cations to models for which the Marshall-Peierls sign rule holds, the theorem 
is exactly obeyed at virtually all levels of implementation in different CCM 
approximation schemes. 

Finally, we are also interested in examining the phase transition points, as 
the strength of the next-nearest-neighbour interactions is varied, at which the 
Neel antiferromagnetic long-range order present in the 2D square-lattice case 
at the Heisenberg point vanishes; and in studying whether there is any rela- 
tionship between the phase boundaries and the points at which the Marshall- 
Peierls sign rule breaks down. We note that any reliable information gained 
on the pattern of the signs of the multi-spin cluster coefficients in the de- 
composition of the ground-state wave function should be very useful for two 
distinct reasons, namely (i) for use in devising improved trial starting wave 
functions for future QMC calculations, and (ii) for spotting possible patterns 
for the signs of the cluster coefficients in different phases or different regimes 
of coupling constants. The latter could be used, inter alia, to suggest possible 
generalisations of the Marshall-Peierls sign rule, and thereby to motivate the 
search for the proofs of suitably generalised theorems. Any such generali- 
sations would clearly have immediate impact for a next generation of QMC 
calculations. 

The rest of this paper is organised as follows. In Sec. II we discuss the 
Marshall-Peierls sign rule and the CCM formalism. The sign rule is first 
outlined for the Heisenberg model on bipartite lattices, before we describe 
the J1-J2 model and previous results. The CCM formalism is then reviewed 
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in very general terms before describing one means of applying it to the J1-J2 
model. Results for the ground-state energy and staggered magnetisation are 
discussed in Sec. Ill, together with results on the breakdown of the Marshall- 
Peierls sign rule for the model. Finally, we present our conclusions in Sec. 



II. THE MARSHALL-PEIERLS SIGN RULE AND THE CCM FORMALISM 

A. The Sign Rule For The Spin-Half Heisenberg Antiferromagnet 

In this section we consider the spin-half Heisenberg antiferromagnet 
(HAF) on a bipartite lattice, where the Hamiltonian is given by, 
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and the index i runs over all lattice points and j runs over the nearest 
neighbours to i. The angular brackets indicate that we count each nearest- 
neighbour bond only once. We note that for a bipartite lattice we can divide 
the lattice into two sublattices such that if i is on one particular sublattice 
then j must be on the other and vice versa. 

For the one-dimensional (ID) linear chain, there is an exact solution to 
this model via the Bethe Ansatz technique [1]. For the two-dimensional (2D) 
square-lattice HAF, there is no exact solution to this problem, though many 
approximate calculations, including those using various quantum Monte Carlo 
[3,4] (QMC) methods and exact series expansion [34] techniques, have been 
performed. 

Although no exact solution is known for the 2D case stated here we note 
that there is an exact sign rule first derived by Marshall [5] (and which we 
shall refer to here as the Marshall- Peierls sign rule). The rule for the square 
lattice HAF is in fact an illustration of the more general Marshall-Peierls 
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sign rule for the HAF on any bipartite lattice. This sign rule provides exact 
information regarding the signs of the expansion coefficients of the ground- 
state wave function in an Ising basis, which is denoted {|/)}. The exact 
ground-state wave function for an AT-body spin system may be written as, 

i*)=5:*,i/), 

where are the expansion coefficients. We now divide the bipartite lattice 
into its two sublattices, denoted A and B, such that each nearest-neighbour 
site to an A sublattice site is on the B sublattice and vice versa. If the number 
of up spins on the A sublattice is called pi then it is possible to show [5] that 
the coefficients satisfy 

*7 = i-lf'ai, 

where the new coefficients aj are all positive. This exact information regarding 
the signs of the coefficients may be used to define the nodal surface of the 
wave function in this basis, and hence is of use in QMC calculations [3,4]. 

B. The J1-J2 Model 

We shall now discuss the spin-half J1-J2 model on the ID linear chain 
and the 2D square lattice. The Hamiltonian is given by, 

N N 

H = Si.Sj + J2 ^ Sj.Sfe , 

where the sum on {i,j) runs over all nearest-neighbour pairs of sites, counting 
each pair (or bond) once and once only; and the sum on ((i, k)) similarly runs 
over all next-nearest-neighbour pairs of sites, again counting each pair (or 
bond) once and once only. We note that in order to consider a wide range 
of the coupling parameters Ji and J2, it is useful to introduce the variable lo 
such that Ji = cos uj and J2 = sin uj. 
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In ID, no exact solution has been found for general values of the coupling 
constants J\ and J2, though there are some exact solutions including the 
Heisenberg point (J2 = 0) and a point at J2/ J\ = 0.5 at which the ground- 
state is fully dimerised [35]. Previous coupled cluster method (CCM) [31-33] 
and density matrix renormalisation group (DMRG) [9] calculations have very 
successfully been carried out for this model. The phase diagram is compli- 
cated, with three distinct phases. These phases may be characterised for our 
purposes as ferromagnetic, antiferromagnetic, and frustrated. The ferromag- 
netic phase is a highly degenerate phase in which the ground-state energy is 
equal to that of the classical fully aligned state. There is a first-order phase 
transition at Ji = with negative J2 to an antiferromagnetic phase. The an- 
tiferromagnetic phase classically has its energy minimised by the Neel state, 
and the quantum-mechanical phase transition point to the frustrated phase 
is at (or is very near to) J2/J1 = 0.5. The frustrated phase classically con- 
tains a spin 'spiral' state which has a periodicity which varies with the ratio 
of the coupling constants J^jJx- There is some evidence that this changing 
periodicity with J2I J\ might also be seen in the quantum- mechanical system 
[33,36]. 

For the square lattice there are no exact results, though approximate 
spin wave theory (SWT) [37] calculations, exact diagonalisations of finite- 
sized lattices [38], and CCM [31] calculations have been performed. The 
ferromagnetic to antiferromagnetic phase transition point is, as for the ID 
case, at Ji = with negative J2, and the antiferromagnetic to frustrated 
phase transition is believed to be near to J2/J1 = 0.5. 

The Marshall-Peierls sign rule, as discussed in Sec. II.A, is true for the 
Heisenberg model on a bipartite lattice. It is simple to prove that it is also 
preserved for the J1-J2 model with negative J2 and positive J\. However, it 
is not in general true for positive J2 and positive J\. In fact, the results from 
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ID short-chain calculations [36] suggest that the breakdown occurs very near 
to the Heisenberg point, at J2/J1 = 0.032 it 0.003. By contrast, finite-size 
lattice calculations [38] on the square lattice indicate that the sign rule at the 
Heisenberg point may well be preserved up to some critical value of J2I J\ in 
the 0.2 < J2/J1 < 0.3. 

C. The CCM Formalism 

In this article we wish to perform CCM calculations for the J1-J2 model 
in the antiferromagnetic regime. We now present a brief survey of the CCM 
formalism and note that a much fuller account of the formalism as applied to 
quantum spin-lattice problems has been given in Ref. [8]. A more extensive 
overview of the method and its applications has also been given in Ref. [17]. 
The starting point for any CCM calculation is the choice of a normalised 
model or reference state, denoted |^>). We define a complete set of mutually 
commuting, multi-spin creation operators {C/} with respect to |^>) such that 
any Ising state [/) may be obtained as 

|/)^C+|<I>), 

for an A^-body spin system. The ground-state wave function has previously 
been written in Eq. (||) as a linear combination of the states {[/)}, and we 
now introduce the usual CCM parametrisations of the ket and bra ground 
states which are given by. 

The ket-state correlation operator in Eq. (P) is, as we can see, formed from a 
linear combination of the creation operators {C/} multiplied with the relevant 
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ket-state correlation coefficients {Sj}. The Hermitian adjoints of the multi- 
spin operators {Cf} are the multi-spin destruction operators {C^}, and the 
bra state in Eq. (^) is formed by the linear combination of these destruction 
operators multiplied with the corresponding bra-state correlation coefficients 
{Sj}. The bra and ket states, defined by Eqs. (|6|) and (^), are not manifestly 
Hermitian adjoints of each other and so the variational property of an upper 
bound on the ground-state energy is not preserved. However, we note that the 
Hellmann-Feynman theorem is preserved. We also note that since {^\Cf = 
= Cyl^) by definition, we have the explicit normalisation relations, (^|^') = 
(§1^) = ($1$) = 1. 

The ground-state expectation value of the energy may now simply be 
written using the Schrodinger equation, H\^) = Eg\^), as. 

Eg = {^\e-^He^\^) . 

Equation (|8|) shows an example of the well-known similarity transform which 
plays a crucial role in the CCM formalism. We further note that the similarity 
transform of any quantum mechanical operator may be written in terms of a 
series of nested commutators, so that for the Hamiltonian H we have 

e-'He' = H+[H,S] + ^[[H,S],S] + --- . 

The infinite series of Eq. (^) terminates at finite order if the Hamiltonian 
H contains sums of products of only finite numbers of single-body operators, 
as is almost always the case and is, indeed, true for the model considered 
here. We also note that each time we perform a commutation operation 
in Eq. we produce a link or contraction, so that every single operator 
in each S within the nested commutator expansion is directly linked to an 
operator in the original Hamiltonian. In this way the Goldstone linked cluster 
theorem is satisfied and the expectation value of the energy, as well as all 



other expectation values, are size-extensive (i.e., they are well defined in the 
asymptotic thermodynamic limit ^ oo at all levels of approximation for the 
operator S). Indeed, the CCM works from the outset in the thermodynamic 
limit. 

We now wish to find values for the ket-state and bra-state correlation 
coefficients. We do this by defining the expectation value, H = {^\H\^), and 
by requiring that this quantity is a minimum with respect to the ket-state 
and bra-state correlation coefficients. Hence, we have, 

6H/6S1 = ^ ($|(77e-^i7e^|$) = , V/ ; 
6H/dSi = ^ {<^\Se-^[H, C/]e^|«>) = , V/ / . 

This formalism is exact in the limit that we include all possible multi- 
spin cluster correlations within S and S, though in any real application this 
is usually impossible. We therefore need to consider approximation schemes 
whereby the expansions of S and S in Eqs. (^) and (^) may be truncated 
to some finite or infinite subset of the full set of independent (fundamental) 
multi-spin configurations. The three most commonly employed schemes have 
been: (1) the SUBn scheme, in which all correlations involving only n or fewer 
spins are retained, but no further restriction is made concerning their spatial 
separation on the lattice; (2) the SUBn-m sub-approximation, in which all 
SUBn correlations spanning a range of no more than ni adjacent lattice sites 
are retained; and (3) the localised LSUBm scheme, which retains all multi- 
spin correlations over distinct locales on the lattice defined by m or fewer 
contiguous sites. 

In the next subsection we consider the application of the CCM to the 
J1-J2 model in the antiferromagnetic regime. 



12 



D. The CCM Applied to the J1-J2 Model 



As stated in the previous section, the starting point for any CCM calcu- 
lation is the choice of the model (or reference) state. Here, we choose the 
classical Neel state to be our model state, in accordance with previous CCM 
calculations [8,31], in order to study the antiferromagnetic regime of the Ji- 
J2 model. We visualise the Neel state by again dividing the lattice into two 
sublattices, A and B, on which each of the nearest-neighbours sites to a given 
sublattice site are on the other sublattice. We populate the A sublattice 
with 'up' spins (i.e., eigenvectors of the operator with eigenvalue +^) and 
the B sublattice with 'down' spins (i.e., eigenvectors of the operator with 
eigenvalue — ^). 

In order to perform a CCM calculation we would like to treat each site 
on an equal footing. We do this by performing a rotation [8,31,33] of the 
local axes of the spins on the A sublattice ('up' spins) by 180° about the 
y-axis such that all spins on each sublattice appear mathematically to point 
downwards (i.e., in these new local axes). Since this rotational transformation 
is a canonical one, it has no effect on the commutation relations. It does 
however have a number of consequences. Firstly, the Hamiltonian is re-written 
in the local coordinates as. 



h=-j,y: 



+J2 E 

{{i,k)) 



We also note that the set of creation operators {Cf} may now be formed 
purely from products of spin raising operators with respect to the rotated, 
'ferromagnetic' model state. We write this expression for an Z-spin cluster 
as = sf^s^^...sf^. Conversely, the destruction operators are now formed 
purely from the spin lowering operators in an analogous manner, where CJ" = 
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The Marshall-Peierls sign rule for the Heisenberg model is also modified. 
We obtain a new and exact rule for the Hamiltonian of Eq. (|l2|) in an ex- 
pansion of the ground-state wave function in terms of an Ising basis, {|/)}, in 
the local, rotated spin coordinates. The corresponding expansion coefficients, 
{^'/}, must now be positive for all of the states labelled by I. (A proof of 
this statement is not given here, but it is made in exactly the same manner as 
that of Marshall [5].) The {^i} coefficients are henceforth explicitly stated 
in relation to the Ising basis in the local, rotated spin coordinates. 

We now wish to provide a link between the {^i} coefficients, in terms of 
the local axes, and the CCM ground-state parametrisation of the ket state 
of Eq. (^). This is done by applying the destruction operator , for a 
particular cluster defined by the index /, to the expressions for the ket-state 
wave function of Eqs. (^) and (^). Note we choose only one ordering out of 
the indices {ii,i2, .-.,11} of the total of N{l\)i> possible equivalent orderings 
for CY on the lattice, where is a symmetry factor dependent on the lattice. 
We therefore write the {^/} coefficients as, 

= (<I>|C7e^|<I>) ^ ms~s-...s- e'm . 

Note that Eq. (^) contains the implicit assumption that the spin raising op- 
erators in of Eq. @, which are used to define |/) with respect to l^*), have 
only one ordering with respect to permutations of the indices {11,12, 

Again, it should be noted that in practice one restricts the choice of the 
clusters contained within S to some well-defined approximation scheme. To 
keep the calculations as self-consistent as possible, we restrict the choice of 
the {^'/} coefficients to be for only those Ising states defined in Eq. (^) which 
correspond to the clusters used in S. 

In the next section we describe our results for the ground-state expectation 
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values for high-order, approximate CCM calculations which are determined 
computationally [8]. We also detect critical points in the CCM equations 
which are taken to be signatures of phases transitions in the real system. 
Once the ket-state correlation coefficients are found it is then possible to 
obtain approximate results for the {^'/} coefficients, again via a computa- 
tional approach, and we discuss CCM results concerning the breakdown of 
the Marshall-Peierls sign rule function of J2/ Ji • 

III. RESULTS 
A. Ground-State Expectation Values 

The ground-state energy of Eq. ^ is approximately obtained once the 
CCM equations are first derived and then solved for a particular approxima- 
tion scheme and approximation level. Descriptions of the method are given 
in Refs. [31-33]. Details of how one may obtain a computational solution for 
high-order LSUBm approximations is given in Ref. [8]|^. We simply quote 
the results here for this model using the Neel model state and the interested 
reader is referred to these articles. 

The LSUBm results for the ground-state energy of the ID J1-J2 model 



^It should be noted that the calculation of Ref. [31] was mostly concerned with SUB2 calculations 
for the spin-half J1-J2 model. However, a calculation for the square lattice system in which only 
nearest-neighbour correlations and four-body correlations between four spins on the unit square 
were retained was also performed. This calculation was referred to as an 'LSUB4' calculation within 
the text in this reference to denote the addition of the extra, single type of four-body correlation. 
However, this 'LSUB4' calculation was not the same as the LSUB4 calculation which we perform 
here which now contains all two-body and four-body correlations in a locale defined by m=4:. 
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converge very well over the range —tt/2 < lo < tan~"^(0.5). We note that the 
LSUBIO results agree to within 1% of those obtained by extrapolating the 
results from exact diagonalisations for short chains [31] over this range, though 
we do not provide a plot of this here. In 2D, we see in Fig. Q that our results 
are again extremely well converged over the range — 7r/2 < u < tan^-'^(0.5). 
In Table | results are given for the ground-state energy of the square lattice 
system as a function of ^ = tanw for —0.5 < < 0.5 for the LSUB6 and 
LSUB8 levels of approximation. 

We note that in 2D the CCM results for the ground-state energy display 
characteristic terminating points at certain critical values of uj. At these 
points the second derivative of the ground-state energy with respect to to 
may also be determined, and we note that at these critical values of uj this 
quantity diverges. This type of behaviour has been observed previously [8] 
and is associated with a phase transition in the real system. The critical 
value in 2D near to the ferromagnetic phase transition, denoted ujp, is given 
in Table We see that the LSUBm results are clearly converging to the 
exact value of lof = — vr/2. It is known [31] that the SUB2 approximation 
predicts this point exactly in both ID and 2D. 

As is seen from the entries in Table II, the antiferromagnetic point, de- 
noted uja, is detected in 2D with the LSUB6, LSUB8, and SUB2 approxima- 
tions. It is not observed within the LSUB4 approximation. We can see that 
the LSUBm critical value of lo decreases with increasing truncation index m, 
and an simple extrapolation [8] in the limit m ^ oo gives a value for the 
phase transition point of ^ ~ 0.61. 

We now introduce the sublattice magnetisation, which characterises the 
degree of quantum order inherent in the CCM wave functions. By inserting 
the CCM parametrisations of Eqs. (^) and we find. 
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r) N ^ N 

M^--Y.{^\sl\^) = --Y.{'f\Se"'sle'\'^) , (14) 

k=i k=i 

where s| is in the local coordinates of each sublattice. Evaluation of the sub- 
lattice magnetisation requires both the ket- and bra-state cluster correlation 
coefficients. The actual procedure to do this is straightforward, and is also 
described in more detail elsewhere [8]. 

The sublattice magnetisation in ID is non-zero in the range — 7r/2 < uj < 
tan^"'^(0.5), though we note (see [31]) that it is greater than zero but mono- 
tonically decreases with increasing LSUBm approximation level for all to in 
this range. The sublattice magnetisation is zero in the true solution of this 
model in ID, and although the CCM LSUBm (and SUB2 [31]) results are 
non-zero we expect that with increasing level of LSUBm approximation this 
would be better reflected in the CCM solution. 

Figure |2| illustrates that that the situation is much clearer in 2D. We can 
see that the sublattice magnetisation in Fig. |2| is converging to a non-zero 
value over essentially all of the range — 7r/2 < lo < tan^^(0.5). We note that 
there are divergences in the sublattice magnetisation which are observed at 
precisely the same points as the critical points, up and uja, of the energy in 
2D. This reinforces our conjecture that these critical points are reflections of 
phase transitions in the real system. 

B. The Breakdown of the Marshall-Peierls Sign Rule 

We now consider the Marshall-Peierls sign rule for the spin-half J1-J2 
model on the linear chain and the square lattice. We need to obtain the {^1} 
coefficients, either analytically or computationally, in terms of the ket-state 
coefficients. We then solve the SUB2 or LSUBm equations in order to obtain 
the ket-state correlation coefficients and hence to obtain approximate values 
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for the {^/} coefficients. Note that these calculations are approximate in 
the sense that we only retain certain correlations in S with a well-defined 
approximation scheme, though we are already working in the infinite lattice 
or — > oo limit. 

We note that at each order of LSUBm approximation it is possible to 
perform this process of matching the terms in e'^ in Eq. (|l^ ) to the configu- 
ration CY analytically. In Appendix A we present the exact form of the {^i} 
coefficients in terms of the ket-state coefficients within the LSUB4 approxi- 
mation scheme for the ID linear chain. Furthermore, we see that each of the 
{^i} coefficients corresponding to two-body correlations with respect to [<!>) 
via Eq. (^), is exactly equal to the corresponding CCM two-body ket-state 
correlation operator {Sj} for any approximation scheme used in S. This may 
be proven by noting that there are no one-body spin correlations allowed in 
S (in order to preserve the conserved quantity Sj^ = J2f ^1 ) ) s-i^d then by 
considering the series expansion of the exponential in Eq. (^). 

To obtain higher-order {'I'/} coefficients using Eq. ( [T3| ) we may con- 
veniently use a computational approach. This amounts to partitioning the 
configuration in CJ = s^_^s^^...s^^ into the multiples of S in the series expan- 
sion of e*^ of Eq. (|T3|). It is then possible to identify by simple computer 
algebra the configurations of the partitioned pieces (each referring to an S 
in the series expansion of the exponential), and find a numerical value for 
the {^i} coefficients once the CCM ket-state equations have been solved at 
specific values of w. 

In both ID and 2D, we find that (at all levels of approximation) the 
Marshall-Peierls sign rule is preserved for uj = — 7r/2|+ (i.e., in the antiferro- 
magnetic regime). That is, all of the {^/} coefficients are found to be positive. 
The sign is broken for uj = — 7r/2|_ (i.e., in the ferromagnetic regime) at which 
point at least one of the coefficients becomes negative. Note that the crossover 
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occurs exactly at the phase boundary uo = — 7r/2, and that this is one case 
where the breakdown of the Marshall-Peierls sign rule occurs at exactly the 
same place as the phase boundary. We note that there is a first-order phase 
transition at this point. 

In ID, we find that all the {^'/} coefficients are positive at the Heisenberg 
point for the SUB2 scheme and for LSUBm schemes with m < 8. Above 
this LSUBm level of approximation (i.e., for m > 8) a few of the {^7} coef- 
ficients (which are very small in magnitude) become negative. However, for 
example, at the LSUB12 level we find that these same coefficients, which are 
negative at the LSUBIO level, again become positive, though we find that 
other new coefficients (which are similarly small in magnitude) then become 
negative. These would in turn presumably become positive at a still higher 
level of approximation. Thus, the CCM is completely consistent with the 
sign rule for the Heisenberg model in ID. This is encouraging as this model 
is quite challenging for the CCM with this model state, as our results for 
the sublattice magnetisation have shown. We may also compare the ratios 
of the magnitudes of the {^/} coefficients to those obtained via short-chain 
calculations, as shown in Table (Note that we examine the ratios to elim- 
inate short-chain normalisation factors.) We can see that the correspondence 
between CCM and short-chain calculations is good, though it appears that 
the CCM results are better converged at the LSUBIO and LSUB12 levels of 
approximation than those from the 12-spin and 16-spin chains. Short-chain 
calculations [36] indicate the breakdown of the Marshall-Peierls sign rule at 
J2/ Ji = 0.032 ±0.003, though our CCM results cannot give an accurate value 
for this breakdown point. 

For the square lattice, the situation is found to be much clearer. The signs 
of the {^'/} coefficients are found to be positive at the Heisenberg point at all 
orders of LSUBm approximation and also from the SUB2 approximation. A 
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clear transition from all of the coefficients being positive to one of them becom- 
ing negative is seen, and we believe that this clearly indicates the onset of the 
breakdown of the Marshall- Peierls sign rule. The points at which the LSUBm 
approximation predicts a breakdown, denoted WAf, are shown in Table ||. 
We can see that a simple extrapolation of these points gives a value for the 
breakdown of the sign rule to be at J2I J\ ~ 0.26 which is in good agreement 
with the SUB2 result of J2/J1 = 0.2607 and exact diagonalisations of finite- 
sized lattice calculations [38] which give a corresponding value of J2/J1 in the 
range of 0.2-0.3. At SUB2 and LSUB4 levels of approximation it is the } 
coefficients for the two-spin cluster with separation x + 2y that becomes neg- 
ative first, where x and y are the unit vectors along the perpendicular axes of 
the square lattice. Higher-order CCM LSUBm calculations predict, however, 
that the coefficients for other higher-order, highly-disconnected configurations 
first become negative at a slightly lower value of cj than the coefficient for this 
two-body configuration. This therefore indicates that higher-order multi-spin 
configurations might as be important as this two-body configuration in the 
breakdown of the sign rule for this model. 

We note that the square-lattice results predict that the breakdown of the 
sign rule (wAf ) occurs at a smaller value of at a particular approximation 
level, than the antiferromagnetic critical point [uja) predicted by the CCM. 
In other words, the CCM results predict that there is a region of the an- 
tiferromagnetic regime in which the Marshall-Peierls sign rule is not being 
obeyed. 

IV. CONCLUSIONS 

The CCM applied to the ID spin-half J\-Ji model gives encouraging 
results for the ground-state energy, the critical points and the sublattice mag- 
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netisation. In this paper we have shown how the recent advances in the com- 
putational implementation of the method enable us to obtain useful results 
for the 2D model as well. We find that the 2D case is in many ways simpler 
than the ID case, with more clearly defined critical points and a non-zero 
sublattice magnetisation. 

We have also investigated the relation between the CCM and the Marshall- 
Peierls sign rule. In ID we find the CCM results are consistent with the exact 
results in the region where these apply. In 2D we are able to obtain results 
which are better than the finite-size extrapolations and we can predict the 
point at which the sign rule fails. Our results indicate that this occurs at 
a different point than the phase transition from the simple antiferromagnet 
to the more complicated frustrated phase. We believe that the use of the 
exact Marshall-Peierls sign rules, extended by the CCM method into regions 
where it is not exact, can shed new light on the behaviour of this type of 
quantum spin system. Furthermore, it may provide information about the 
nodal surface that can be used for accurate QMC calculations in a much 
wider range of quantum spin systems than previously. 
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APPENDIX A: ANALYTIC EVALUATION OF THE COEFFICIENTS 



We now present analytical expressions for the {^i} coefficients in terms 
of the ket-state correlation coefficients in S for the LSUB4 approximation in 
ID. The ket-state correlation operator for the LSUB4 approximation in ID is 
given by, 

N N N 

s = biJ2 44+1 + b3j2 44+3 + 54 ^ 44+i4+24+3 > 

i i i 

where the index i runs over all points on the linear chain. To obtain the {^i} 
coefficients we now must choose the configurations in Eq. (|T^): for the 
nearest-neighbour, two-body coefficient, which we shall denote as ^'i, we use 
Ci = sjs~_^i; for the third-nearest-neighbour, two-body coefficient, which we 
shall denote as ^^21 we use C2 = sjsj^^; and for the four-contiguous spin con- 
figuration coefficient, which we shall denote ^'3, we use = sJ sJ_^^^sJ^2''^J+3- 
The result is therefore, 

*2 = 63 ; 

^3 = bl + hbs + 54 ■ 

The values of the coefficients are independent of j, due to the translational 
symmetry of the lattice, and so the index j is chosen arbitrarily from any of its 
N possible values in order to obtain Eqs. (A2-A4). Higher-order LSUBm ap- 
proximations can be handled analogously by making use of computer-algebraic 
techniques. 
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FIGURES 
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FIG. 1. Results for the COM ground-state energy of the spin- half J1-J2 model on the 2D 
square lattice, using the LSUBm approximation based on the Neel model state, with m = 2, 4, 6, 8. 
LSUBm critical points, ujp and ua, are indicated by the boxes. Note that Ji = cos u and J2 = sin u. 
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FIG. 2. Results for the CCM sublattice magnetisation of the spin-half J1-J2 model on the 2D 
square lattice, using the LSUBm approximation based on the Neel model state, with m = 2, 4, 6, 8. 
Note that Ji = cos ui and J2 = sin w. 
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TABLES 



TABLE L Results for the CCM ground-state energy of the spin-half J1-J2 model on the 2D 

square lattice as a function of J2/J1 = tan cu, using the LSUBm approximation based on the Neel 
model state, with m = 6, 8. 



J2/J1 ~ tan(w) 


Eg/N (m=6) 


Eg/N (m=8) 


-0.50 


-0 88237 


-0 88353 


-0.40 


-0.83785 


-0.83902 


-0.30 


-0 79393 


-0 79510 


-0.20 


-0.75071 


-0 75188 


-0.10 


-0.70833 


-0.70951 


0.00 


-0.66700 


-0.66817 


0.10 


-0.62699 


-0.62816 


0.20 


-0.588G8 


-0.58988 


0.30 


-0.55271 


-0.55397 


0.40 


-0.52012 


-0.52164 


0.50 


-0.49311 


-0.49551 
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TABLE II. Results for the CCM critical points of the spin-half J1-J2 model on the 2D square 
lattice. The ferromagnetic and antiferromagnetic critical points are denoted cof and ua respectively. 
The point (lu > 0) at which the CCM predicts that the Marshall-Peierls sign rule breaks down 
for the square lattice is denoted lom- Corresponding results for {J2/Ji)\a and {J2/Ji)\m are also 

sh(nvu, where; J2/./1 = tan 









{J2/Ji)\a 




{J2/Ji)\m 


SUB2 


-7r/2 


0.6416 


0.7470 


0.255 


0.2607 


LSUB4 


-1.702 






0.331 


0.344 


LSUB6 


-1.628 


0.583 


0.660 


0.290 


0.298 


LSUB8 


-1.603 


0.566 


0.636 


0.275 


0.2<S2 


LSUBoo 


-1.572 


0.544 


0.605 


0.255 


0.261 



TABLE III. Results for the ratios of the magnitudes of the CCM LSUBm '^i coefficients, for the 
spin-half, ID Heisenberg model, compared to ratios of the magnitudes of the equivalent Ising expan- 
sion coefficients determined by finite-size, exact diagonalisations. The ^'i coefficient is associated 
with the nearest-neighbour, two-body correlation with respect to the model state |$); and ^^2, *3, 
and ^'4 arc, respectively, the corresponding coefficients associated with the 3'^'^-nearest-neighbour, 
two-body configuration, the four-contiguous-spin configuration, and the six-contiguous-spin con- 
figuration. 



Ratio 


12 Spins 


16 Spins 


LSUBIO 


LSUB12 


l*3l 
l*l| 


0.7436 


0.7163 


0.6720 


0.6758 


l*2| 
l*l| 


0.1674 


0.1552 


0.1413 


0.1381 


l*4| 


0.4850 


0.6155 


0.5157 


0.5248 
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